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Abstract 

We study semileptonic B decay to the exclusive charmless states it, p/u, n and rf using the 
16 fb- 1 CLEO T(45) data sample. We find B(B° -► 7r~£ + u) = (1.37 ± 0.15 sta t ± 0.11 sys ) x 10~ 4 
and B(B° -»• p~l + v) = (2.93 ± 0.37 sta t ± 0.37 sys ) x 10" 4 , and find evidence for B + -► n'i + v, with 
B(B+ -► r/'£+i/) = (2.66±0.80 sta t±0.56 sys ) x 10~ 4 . From our B -► vrfz/ rate for g 2 > 16 GeV 2 and 
lattice QCD we find \V ub \ = (3.6 ± 0.4 stat ± 0.2 syst +°;° t ) x lO" 3 . 
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The magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [T] element V u b plays 
a central role in tests of the matrix's unitarity, which is the main thrust of the current heavy- 
flavor program. Violation of unitarity would signal existence of new classes of fundamental 
particles or forces. Robust determination of \V u b\ has been the subject of considerable 
theoretical and experimental effort for well over a decade, and remains one of the highest 
priorities of flavor physics, yet it remains one of the most uncertain parameters of the matrix. 
A recent determination, dominated by inclusive measurements, has marginal agreement 
with other inputs to the unitarity constraints [2]. Exclusive charmless semileptonic decays, 
particularly B — > ixlv, provide an alternate route to determination of \V u b\. 

This letter presents a study of the decays B — > ixtv, B — > piu and B — > v^'lv based 
on the full 15.4 million BB data sample collected at the T(4S) with the CLEO II, II. 5, 
and III detectors [31 SI E] at the Cornell Electron Storage Ring (CESR). These results 
include an important crosscheck of the leading measurement of the critical B° — > 7r~£ + z/ 
branching fraction [6] using a different analysis technique in a symmetric beam environment. 
We also present the most precise study of B — > piu, which is an important background 
process to B — ■> ir£v. The results include the first partial rate measurement for a region of 
B — > piv phase space (helicity angle over 90°) that provides about half of the total B — > ptv 
background in the B — > -kCv signal region. This benefits both this and other B — > iriu 
measurements that have a background from B — > plv. Finally, we find first evidence for the 
decay B — > r]lv. This channel is predicted to be enhanced relative to B — > r]lv through 
coupling to the singlet component of the r{ [7J [S] , which could resolve the unexpectedly large 
decay rate for B — > rj'K. These results supersede those of references [HJ [TO] . 

Hadronic form factors (FFs) present challenges experimentally and theoretically for de- 
termination of | V u b | • For the decay B — > V u lv^ where V u is a charmless vector meson, the 
partial width is 



dT 2 
nkq 



S 2 H 2 + °- H + + °+ H - 



dq 2 dC e L~"" U 2 J " (1) 

\V h\ 2 G 2 

Here, k = [ 2 gJ 3M a , k is the V u momentum, Cq (Se) is the cosine (sine) of the angle 9wi 

between the charged lepton in the virtual W (W*) rest frame and the W* in the B rest 
frame, and C± = 1 ± Cg. H± and H are the magnitudes of the W helicity amplitudes, 
which can be expressed in the massless lepton limit in terms of three g 2 -dependent FFs [UJ . 
For a final state pseudoscalar meson P u , H± = 0, and the rate depends on a single FF 



/+ W) ■ (2) 



dT(B^PJu) 2 G 2 F 3 

dq* 1 ubl 2^ 

Experimental measurement of the branching fractions requires knowledge of the q 2 - 
dependence of the FFs. All current measurements follow the previous CLEO analysis [10] 
and mitigate FF shape uncertainties by measuring partial rates in regions of phase space. 
Subsequent determination of \V u b\ requires theoretical input on the FF normalization. Re- 
cent calculation of f+(q 2 ) from lattice QCD using dynamical quarks [12] represent a marked 
theoretical advance. 

The analysis uses the missing four-momentum, p m i SS = pcm — Pvisibie to estimate the 
neutrino four-momentum, p u . All CLEO detector configurations [3jmi5] provide acceptance 
over more than 90% of the full 47r solid angle for charged particles (momentum resolution 
of 0.6% at 2 GeV/c) and for photons (average n° mass resolution of 6 MeV/c 2 ). We utilize 
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global event reconstruction and particle identification algorithms, optimized for neutrino 
reconstruction, outlined in detail elsewhere [T3l [14] . 

Signal leptons are identified with over 90% efficiency in the range 1.0 < pi < 3.0 GeV/c. 
We reconstruct hadronic X u candidates 7r°, 7r ± , 77, p°, p ± , u>, and ?/. We require the recon- 
structed 7T° (77) mass in the 77 mode to be within two standard deviations, about 8 (26) 
MeV/c 2 , of the nominal 7r° (77) mass. These candidates are kinematically fit to improve the 
momentum resolution. Candidate rj (u) decays with 7r + 7r~7r° invariant mass within 10 (30) 
MeV/c 2 of the nominal 77 (uj) mass and p candidates with invariant mass within 285 
MeV/c 2 of the p mass are accepted. The decay of the rj' is reconstructed in both the tttt^ 
and 777r + 7r~ final states, with 77 — > 77. For 77' — > 7T7T7, the invariant mass of the reconstructed 
77' must be within 2.5 (2.75) standard deviations of the nominal rj' mass for q 2 < 10 (q 2 > 10) 
GeV 2 . For 77' — > tj-rtt the quadrature sum of the number of standard deviations from the 77 
and 77' masses must be less than 3.75. 

In events with multiple undetected particles, p m i SS represents p v poorly, causing B — > X c lv 
decays to smear into the signal region of (much rarer) B — > X u tv processes, while smearing 
much of the X u iv out. This mechanism provides our dominant background contribution. 
Therefore, we reject events with multiple identified leptons, which are usually accompanied 
by multiple neutrinos, and events with net charge | AQ\ > 0, an indication of missed particles. 
In the pseudoscalar decay modes, we do achieve some additional precision by independently 
considering the more contaminated \AQ\ = 1 sample. We further require that M^ iss be 
consistent with a massless neutrino within experimental resolution. Because our resolution 
on E miss , cj£; mias (~ 0.2 GeV), is larger than that for p miss (~ 0.1 GeV), the M^ iss resolution 
scales as ~ 2.E miss /<7E miss , so we require |M^ iss /2£ , miss | < 0.5 GeV - a requirement that 
is approximately constant in -E m i SS resolution. The criterion was optimized, using only 
independent Monte Carlo samples, for maximum signal significance. 

Candidate B decays are formed by combining the signal lepton and reconstructed neu- 
trino with X u candidates, taking p v = (|p m i SS | ,p m i SS ) for reconstruction. Energy and mo- 
mentum conservation requirements are cast in terms of AE = Ex u + E# + E v — -Ebcam and 

M Mv = \J i?beam — \Px u + Pi + otp v \^ , which peak at zero and M B for signal decays. Because 
the neutrino resolution dominates the AE resolution, we improve the resolution by 

choosing a such that Ex u + Ei + aE v — -Eboam = 0. Before computing the kinematic variables 
q 2 and cos^vy^, the reconstructed neutrino is rotated to force Mmv to the value of Mb to 
optimize resolution in these variables. This procedure produces resolutions of 0.3 GeV 2 and 
0.03 for q 2 and cos^vy^, respectively. 

Continuum e + e~ — > qq, where q = u, d, s, or c, can be distinguished from BB events 
by examining a variety of event shape variables. Specifically, we examine the ratio of the 
second to zeroth Fox- Wolfram moments [T3], the angle between the candidate thrust axis 
and the thrust axis of the rest of the event, the angle between the event thrust axis and 
the beam, and the momentum flow into nine double cones whose axes are centered on the 
beam. A Fisher discriminant [16] is computed using these twelve variables. Requirements on 
this discriminant are tuned for each decay mode and each q 2 and cos Owe region to maximize 
signal significance. Due to their jet-like nature, continuum backgrounds typically reconstruct 
with low q 2 and can be largely isolated at q 2 < 2 GeV 2 . 

We fit seven coarse bins [13] of the AE versus M Mu distributions, separately reconstructed 
within each phase space region summarized in Table |T] to allow extraction of a partial branch- 
ing fraction for each region. The bins span 5.1750 GeV/c 2 < M^ v < 5.2875 GeV/c 2 and 
—0.75 GeV < AE < 0.25 GeV. We are primarily sensitive to signal within the bin defined 
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by 5.2650 GeV/c 2 < M Mu < 5.2875 GeV/c 2 and -0.15 GeV < AE < 0.25 GeV. We also 
use coarse lineshape information to help separate signals involving resonances in the final 
state from background. For the plv (ulv) mode, we accomplish this by separation of the re- 
constructed distributions into three 195 MeV/c 2 (20 MeV/c 2 ) intervals in reconstructed tttt 
(Sir) mass centered on the nominal resonance mass. In 77' — > 7r + 7r~7, we use four M W7r bins 
which span the range from 300 MeV/c 2 to 900 MeV/c 2 . The fit utilizes a binned maximum 
likelihood approach extended to include the finite statistics of the fit components p2] . 

A GEANT-based Monte Carlo [18] is used to model the distributions for the signal, 
b — ► c background, and backgrounds from other B — ► X u £z/ decays that are not explicitly 
fit, e.g. B — > a £zA The signal Monte Carlo (MC) is divided at the generator level into 
the phase space regions of Table |TJ We obtain a full set of reconstructed distributions from 
each which are normalized in the fit to a rate parameter for that region. We fit all modes 
(and their subregions) simultaneously. The data in one phase space region thereby controls 
the normalization for that region's cross-feed into all other regions and modes, minimizing 
our dependence on a priori form factor and branching fraction information. The input 
form factors for the signal MC samples are derived from an unquenched lattice (LQCD) 
calculation [12] for nlv and a light cone sum rules (LCSR) calculation [19] for pjiolv. Data 
taken below the T(4S), smoothed to reduce statistical fluctuations [13], is used to model 
residual continuum backgrounds. Small contributions from fake signal leptons are modeled 
using purely hadronic data and measured lepton fake rates. 

To increase sensitivity to rflv decay, the candidates are divided into two bins in q 2 : one 
greater than 10 GeV 2 , which typically has higher B — > X c lv background, and one less 
than 10 GeV 2 . This division enhances the the significance of the expected signal-rich region 
without introducing form factor dependence associated with a strict q 2 requirement. We fit 
only for the total B — > rjiv rate by using the ISGW2 model [20] to fix the ratio of efficiency 
corrected yields in these two regions. 

We minimize multiple candidates per event to simplify statistical interpretation. Within 
each 7r, 77 or rf — ► r/7r + 7r~ mode, we allow only one candidate within the fit's AE — M^ v 
region. Within a mode, multiple candidates in that region are resolved by choosing that with 
the smallest |Ai?|. For the vector and the 77' — > 7r + 7r~7 modes, the \ AE\ criterion combined 
with large combinatorics outside the signal region induces efficiency loss. To mitigate loss, 
we select the best candidate separately within each M nn or M 3n range. This procedure 
induces only very slight peaking in backgrounds, which is modeled by the MC generated 
distributions used to fit the data. While a given event can still contribute to more than one 
mode (or mass range), tests of the fitter using numerous MC samples as toy data sets verify 
that signal yields and statistical errors are extracted without bias. 

Signal selection efficiencies, averaged over phase space, for the restricted signal region 
range from rs 4% for the 7r iv mode down to ~ 0.5% for the rj'lv mode. 

The nominal fit includes one free parameter for the yield in each of the 4 (5) regions in 
nlv (plv) phase space and free parameters for the total rjlv and rj'lv yields. Within each 
phase space region, the charged and neutral tc and p rates were fixed to be consistent with 
expectations from isospin, B lifetime differences, and relative charged/neutral B production 
at the T(4S) [21]. The total B + — ► uj£ + v rate is constrained to the B + — > p°i + v rate. To 
minimize systematic biases, the normalization of b —>■ c background floats freely for each 
mode and for each | AQ\ sample. The sum of the signal rates and the feed-down from higher 
mass B X u lv rates is constrained by recent inclusive b — > u measurements at the lepton 
endpoint [22] • There are a total of 26 parameters. 
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M h|v (GeV/c 2 ) 

FIG. 1: Mmu for candidates in the AE signal region [-0.15 GeV,0.25 GeV] for the B -> ixlv 
modes combined and summed over q 2 . Dotted lines indicate the binning used by the fitter. Fit 
components, from bottom to top are, b —* c, continuum, fake signal leptons, other B — > X u £z/, 
£? — ► cross-feed, £? — > p^i/ cross-feed, B — > 7r£^ cross-feed, and signal. 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
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FIG. 2: M^,,- in the restricted signal region for the combined B — » modes, with components as 
described in Figure [T] Arrows indicate the range included in the fit. 



Table [I] summarizes our branching fraction results. The fit yielded — 21n£ = 541, and we 
note that the statistical errors on some of the 532 bins are not Gaussian. Figure [T] illustrates 
a Mhe u projection of the data and fit components for the B — ► 7t£u mode. The signal, 
peaking at Mmv = Mb, is clearly visible. We also show M n7T (Figure [2]) for B — > plv. In all 
cases the data are nicely modeled by the fit. 

A detailed summary of the systematic errors can be found in [13] . The dominant experi- 
mental systematic error is neutrino reconstruction efficiency. This uncertainty arises mainly 
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from systematic effects within the MC simulation, and we have estimated it by randomly 
discarding or smearing the reconstructed tracks and showers based on independent studies 
of tracking and calorimeter resolution and efficiency We also vary the B — > X c iv form 
factors and branching fractions within present experimental limits. To estimate systematic 
bias due to non-resonant B — > -nulv contributions, we study reconstructed B — > 7r°7r°£z/ 
candidates. Non-resonant mriu components generated with p lineshapes are included in 
the fit with relative strengths of the various charge combinations constrained by isospin 
and angular momentum considerations [§J. The 7r°7r°^z/ mode then limits the level that the 
p lineshape can be projected out of the (unknown) non-resonant shape. No statistically 
significant non-resonant bias of the B — > piv rate is observed. 

While our fitting approach significantly reduces dependence on a priori knowledge of 
the FFs, residual FF dependence due to efficiency variation within each phase-space region 
remains. We study this dependence with several other FF calculations for B — ► irlu [201 
|23j El] and B — > piv j20j |25l [26] , chosen to span conservative ranges in shape. We vary the 
B — > plv FF keeping the B — > -ntv FF fixed to its nominal, and vice versa. The systematic 
uncertainty is estimated be one-half of the total spread of the results. 

We find evidence for B (B + r]'£ + u) = (2.66 ± 0.80 ± 0.56) x 10" 4 at 3a significance 
and set an upper limit at the 90% confidence level on the branching fraction for the decay 
B (B + -> r]£ + u) < 1.01 x lO" 4 . The latter results require B (B+ -> r]'i + u) /B (B + -> r]£ + u) > 
2.5 at 90% confidence level, an indication of enhanced coupling to the singlet component 
of the rf . From a model dependent study (see [T5]) of the rf^ 1 rates following Beneke and 
Neubert [S], we estimate B{B~ — > K~rf) = (84jl|| x 10 -6 , in good agreement with the 
current experimental value of B(B~ — > K~rf) = (71 ±4) x 10~ 6 [27J. The 3a significance of 
our B (B + — > rj't^v) result is estimated by using a toy MC procedure to combine systematic 
and statistical errors and determine the probability that zero rate would produce a result 
at or above our measured value [T3] . 

Our B — > piv results currently provide the most precise test of the B — > p form factor 
calculations to date. We fit the predicted form factor shape to the measured rate in the five 



TABLE I: Summary of the phase space subregions and the partial and total branching fraction 
results. 
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phase space regions. The resulting fit has four degrees of freedom, and we find the x 2 to be 
4.5, 9.0, and 4.3, for the LCSR [[9], ISGW2 [20], and UKQCD [26] calculations. The LCSR 
and LQCD calculations, extrapolated over all phase space, agree with the data much better 
than the ISGW2 model. 

The full potential for extraction of \V u t,\ will only be realized once techniques that can 
utilize a broad range of q 2 have matured. For now, we extract a value for \V u b\ by combining 
our measured rate for B° — ► 7r~£ + i/ in the q 2 > 16 GeV 2 region with the B° lifetime and the 
partial width prediction, r/|^| 2 , from LQCD [12], which calculates FFs using a minimal 
number of assumptions. We find \V u b\ = (3.6 ± 0.4 ± O^Iq^ x 10~ 3 . The errors are, in 
order, statistical, experimental systematic, and theoretical. 

In summary we have measured the branching fractions B (B° — > tt~£ + v) = (1.37 ± 0.15 ± 
0.11) x 10~ 4 , B (B° -> p~£+v) = (2.93 ± 0.37 ± 0.37) x 10~ 4 , and B (B+ -> r)'£ + v) = (2.66 ± 
0.80±0.56) x 10" 4 . We find \V ub \ = (3.6 ± 0.4 ± 0.2±g;|) x 10~ 3 from B (B° -> ir~£ + v) in the 
q 2 > 16 GeV 2 region, where our results are consistent with other measurements [SI |28j |29] . 
These results constitute one of the two most precise measurements of B (B° — > it~£ + v) [S], 
the most precise measurement of B (B° — > p~£ + v), and the first evidence for B — > rj £v. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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